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Abstract 

The behaviour of a system of two two-level atoms, both identical 
Q_t as well as non-identical, under the influence of a continuously varying 

dipole-dipole coupling parameter for a wide range of other parame- 
ters such as the atomic level separations and the coupling strength 
of the external radiation field is explored. A detailed analysis of the 
behaviour of various level populations and some of the atomic coher- 
ences is presented. The influence of the competing effects between 
the Rabi field strength and the dipole dipole interaction on the two 
photon absorption probability is explored in detail both for identical 
as well as non-identical atoms. There are significant variations in the 
behaviour of identical and non-identical atoms as a function of the 
dipole coupling strength, even for small amounts of non-identity of 
the atoms, which in this study is incorporated through the atomic 
level separations. 

1 INTRODUCTION 

Cooperative effects pQ, for example collective emission, are a general charac- 
teristic of macroscopic systems consisting of a large number of atoms confined 
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in a region of space much smaller than an optical wavelength. However, in 
order to understand the collective behaviour of macroscopic systems, it is 
essential that we have a microscopic formulation of the interaction between 
the atoms and the electrons. Mention has to be made that cooperative effects 
are not important in most experiments involving interaction of atoms with 
laser fields, which are in the dilute gas regime, where the confinement region 
is much larger than an optical wavelength. However, situations such as in 
a Rydberg gas where the dipole moments are very large or in Bose Einstein 
Condensates (BEC) where the densities are higher, the collective behaviour 
of atoms indeed has strong influence on the outcome of the experiment. 

We study a prototype model consisting of a pair of interacting atoms in 
the presence of a radiation field of arbitrary strength. Even though such a 
model may seem rather elementary, the study is significant since this model 
can serve as a basic building block in understanding of the cooperative ef- 
fects among interacting atoms. Recent advances have enabled trapping and 
studying of very few atoms in a trap, in light of which this study is all the 
more relevant. 

A large amount of work exists that study cooperative effects among few 

atoms [TJ[2j[3j|lJ[5l[6j[7J[8j[9j ^ ' °^ P ar ti ciuar interest being the case of a 
single pair of two-level atoms. The justification for the latter comes from the 
fact that collective phenomena basically are a result of pair-wise interaction 
among atoms and hence the study of two atoms forms the basic building block 
in treating the collective interactions. One such earlier work [TT] reports that 
correlation between two identical atoms can give rise to new resonances in the 
resultant fluorescence spectrum. In another study[12j it has been shown that 
dipole-dipole interaction between atoms may lead to the suppression of the 
fluorescence whenever the atoms are placed within one optical wavelength. 
In reference [T3], Kim et al. have shown that cooperative effects can cause 
modification in transition rates and level shifts, thus giving rise to significant 
changes in the multiple jump dynamics. 

In yet another study, Varada and Agarwal [14J have observed an en- 
hancement of the two photon resonance due to the dipole-dipole interaction 
strengths. Another aspect, namely, the entanglement characteristics of such 
systems were explored via a study of two photon correlations by von Zanthier 
et al.pj)]. 

In general, in a collection of atoms, it is observed that, the first excited 
atom prevents the excitation of neighbouring atoms due to the dipole-dipole 
interaction, which is referred to as the dipole blockade [TBI ITT] . This dipole 
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blockade can give rise to the production of single atom excited collective 
states. 

In a more recent work [TS] the effect of dipole blockade on the dynamical 
evolution of two two-level atom system in the presence of a driving field 
was studied to understand the entanglement characteristics of the system. 
Rydberg atoms, by virtue of their large dipole moments exhibit strong dipole- 
dipole interactions which can give rise to entanglement [T^J EDI EI]- This 
interesting property of Rydberg atoms has led to proposals for using them 
in quantum computation schemes |22j. 

In view of the rich and varied applications the system of two two-level 
atoms has to offer, it becomes an interesting subject of study. One ought to 
mention that most of the earlier work reported is carried out in the dressed 
state formalism. In this paper, we examine in detail, the behaviour of the 
two two-level atom system both for identical as well as non-identical atoms. 
We work in the product space of the single atom states as we believe that 
scaling up to more atoms is straight forward. In this basis, the two two-level 
atom system can be modeled as an effective four level atom system. 

The study carried out here is valid for arbitrary strengths of various in- 
teractions namely the electric dipole interaction between the atom and the 
external radiation field, the dipole-dipole interaction between atoms. The 
master equation approach [23J is used to incorporate the spontaneous emis- 
sion decays. 

Density matrix formalism is employed to obtain the equation of motion 
for the various atomic level populations and coherences. The resulting sixteen 
coupled differential equations governing the dynamics of the system can be 
solved both in steady state regime as well as in time domain. Emphasis 
is laid on exploring the influence of competing behaviour between the Rabi 
field strength and the dipole interaction on the level populations, atomic 
coherences and the total one photon as well as two photon absorption cross 
sections, across different frequency regimes. 

As mentioned earlier, Varada and Agarwal[14j have reported the enhance- 
ment of two-photon resonance due to the presence of dipole-dipole coupling 
between non-identical atoms. In this paper, we have extended their study 
to a wider parameter regime to explore the behaviour of this two-photon 
resonance in greater detail. We have observed that the enhancement of the 
two-photon resonance is not a monotonic feature. Instead, beyond a certain 
critical value of the dipole coupling strength, it begins to decrease again. It 
is this feature of the two-photon resonance that has been addressed in detail 
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in this paper. 

The organization of this paper is as follows. In Sec. II, we present the 
mathematical formulation of the problem. The basic equations that gov- 
ern the dynamics of the system of two two-level atoms interacting with the 
single-mode laser field are written down in the rotating wave approxima- 
tion. In section III, numerical results in steady state of various quantities 
are presented for a wide choice of parameters, both for identical as well as 
non-identical atoms. It is seen that a small deviation from identity of the 
atoms gives rise to significant changes in the behaviour of the system. 



2 THE MODEL 

We consider two isolated two level atoms fixed in position with a separation r, 
driven by a single mode resonant laser field, of frequency ojl (wavelength A), 
which is nearly resonant with the transition frequency of one of the atoms, in 
the region r< A. Here |gj) and |e^) are respectively the ground and excited 
states of the i th atom ( i=l,2). The interaction with the field induces dipole 
moment in the atoms, which in turn interact with each other via photon 
exchange. In this model we assume that the external field is propagating 
perpendicular to the inter atomic axis so that the interaction of field with 
the atoms is purely symmetric. We use the semi-classical approach to study 
the problem in which the atomic system is treated quantum mechanically 
while the laser field interaction as well as the dipole - dipole interaction are 
treated classically. As mentioned earlier, the system of two two-level atoms 
can be represented as an effective four level system as shown in figure [TJ 

The product state \g1g2) ((e^)) corresponds to the situation where both 
atoms are in their ground(excited)state. The intermediate state \g1e2) (\e1g2)) 
represents the case when the first (second) atom is in ground (excited) state 
and second(first) atom is in excited (ground) state. For brevity the follow- 
ing notation is introduced, \g1g2) = |l),|ei<7 2 ) = 1 2) , | f?ie 2 ) = |3) and 
\e1e2) = |4). The Hamiltonian for this two two-level atomic system can 
be written as 
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I 9i,g 2 >= |1 > 



Figure 1: Schematic energy level diagram of two two-level atoms as an equiv- 
alent four-level system 

where Wj(i=l,2,3,4) is the frequency of the i th level, ay are the different 
Rabi frequencies which represent different atom-field coupling strengths and 
g represents the strength of the dipole - dipole interaction. 

Several theoretical approaches can be used to treat the system of two 
atoms interacting with a quantized radiation field. One of the more popular 
approaches is the master equation method to obtain the equation of motion 
for the density operator of the atomic system coupled to a reservoir (the 
multimode vacuum field). The master equation approach enables one to 
treat the evolution of the atom plus field system entirely in terms of the 
atomic operators. In this work, we use the master equation formalism to 
study the dynamical evolution of the density operator of the atomic system. 
The equation of motion for the density operator p is given by 

th^ = [H,p]-Cp (2) 

The first term on the right hand side represents the interaction with the 
radiation field and the dipole-dipole interaction whereas the second term 
(The Liouvillean) accounts for different decay mechanisms. Mention has to 
be made that in the equations of motion that are presented in the following, 
the spontaneous emission decays alone have been included. However, the 
equations are suitably modified to take into account the cooperative decay 
wherever in the calculations, the role of the cooperative decay parameter also 
has been examined by suitable modification of the equations of motion. The 
resulting equations of motion for the sixteen density matrix elements in the 
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rotating wave approximation are obtained as follows: 

dtpii = ia* 21 p~2i - ia2iPi2 + i^liPli - icoiph + 2721P22 + 2731/933 
dtP22 = -ia*2iP~2i + ia2iP\2 + ig(ph - P23) + ial 2 Ph ~ i"42P24 + 2742P44 - 2721P22 
dtph = -ialiPii + ia3iP~i3 - ig(ph - P23) + ial 3 ph - ia43P34 + 2743^44 - 2j 31 pl 3 
dtpli = ia 42 p24 + ia43P34 ~ ial 2 ph - ial 3 ph - (2742 + 2743)^44 

d t ph = («Ai - 721)^12 - ia^iPii + ia 2 iP~22 + ia 31 ph - igph ~ ia42Pu (3) 

dtPh = (»A 2 - 73l)pi3 - i a 3lPll + ia 3 lP~33 + ««2lP23 - WP12 - iO-A3P~\A 

dtPu = (iAi + zA 2 - 743 - 742)^14 - ial 2 p~i2 - iahph + ioi 2 \P~2i + ia 3 \P3i 
dtph = (»A 2 - iAi - 721 - 731)^23 + i"2iPi3 - ia^pli + iot* i2 P~i3 - i&A3P~2A + ig(p33 - P22) 
d t p~2A = («A 2 - 721 - 742 - 743)^24 + ia2iP~iA + ia&pii - iot% 2 p~22 - ia^ph + WP~3A 
dtP3A = (iAi - 731 - 742 - 743)^34 + ia3iP~iA + ioi* i3 pii - ial 3 ph - ia^ph + WP~2A 

The remaining equations can be obtained by realizing that 

dtpij = dtp~j*. 

In the above A x = co atom i - u L , A 2 = 0J atom2 - u L , and the total 
detuning A = Ai + A 2 . 

Here 27^ represents the rate of spontaneous decay from level j to i while 
7ij is the decoherence rate. By using the completeness condition pn + p 2 2 + 
P33 + P44 — 1; the 16 equations can be reduced to 15 equations which can be 
cast in the form 

dm 

— = Mm + $ (4) 
at 

where M is a 15 x 15 coefficient matrix and m and $ are column vectors each 
of length 15 which are defined in the following. 

* = [Pll Pl2 Pl3 Pl4 P21 P22 P23 P24 P31 P32 P33 P34 P41 P42 P43] T (5) 

$ = [0 2 742 ia* 42 2 743 ia* 43 - ia 42 - ia 43 } T . (6) 
The steady state solution of the density matrix elements may be obtained 

as 

^ ss = y(t -> 00) = — M _1 $ (7) 

In the next section we present some of our numerical results of various 
level populations and coherences in the steady state, for both identical as 
well as non-identical atoms. With no loss of generality, we assume that 
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the coupling between the atom and the field to be same for both atoms, 
which corresponds to all ck^'s being equal which is denoted by a and all the 
spontaneous decay rates also are taken to be equally = 7) 



3 RESULTS AND DISCUSSION 

The steady state behaviour of level populations, various atomic coherences 
and absorption cross sections has been studied numerically for a wide range 
of parameter regimes, both for identical and non-identical atoms. In this 
section, some representative numerical results for typical values of the pa- 
rameters are being presented. 



3.1 Identical atoms 



The steady state populations pa, « = 1,2,3,4 for identical atoms, are plotted 
as a function of the total detuning A = Aj.+ A 2 , in figure|2]for a given value of 
the Rabi frequency and for different values of the dipole coupling parameter 

P22,P44 and 



2(b), 2(c) and 2(d) respectively show p 



11, 



g. Sub-figures 2(a) 
P33. Each of the sub-figures corresponds to different values of the dipole 
coupling parameter namely g = 0,5, 10 and 20 from top to bottom. As 
already mentioned, the value of Rabi frequency is held constant for all of 
them at a = 2. It is to be noted that all parameters for all the data presented 
in this paper are normalized to 27, the population decay rate. 

When there are no dipole-dipole interactions, i.e., g = 0, the two atoms 
are independent of each other and hence show single resonant peak at A = 
0, characteristic of a single two- level atom (topmost curves in all the four 



subfigures 2(a) 2(d)). 



All possible configurations, such as both atoms being in their respective 
ground state (pn), any one of them in the excited state while the other 
in ground state (p 2 2 and P33) and both atoms in excited state (/O44), are 
all equally probable and hence show a dip/peak value of 0.25 at resonance. 
The resonant dip in p n is matched equally by the increase in p 22 , P33 and 
P44. The rest of the curves in figure^] correspond to non zero values of g in 
increasing order from top to bottom as indicated in the figure. The salient 
features observed in these curves are as follows: (i) The central resonant peak 
decreases in height (depth in case of pn) with increase in g. (ii) A new side 
band appears at a non zero detuning A. Though it is very difficult to solve 
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the equations analytically, we have nevertheless been able to deduce, from 
the numerical values computed over a wide range of parameter values, that 
this peak appears approximately at A = 2a/ (\g\ 2 + 2a 2 ). (iii) This side band 



does not exist for ( 2(c) ), although a small protuberance can be seen for 



intermediate values of g. (iv) The populations P22 and P33 are identical, as 
expected, indicating that single atom excitation does not distinguish between 
the two atoms, (v) With increasing g, the value of the resonant peak/dip at 
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A = identically decreases for all the four states. This re-affirms that (a) 
the energy levels are shifted due to the dipole-dipole interaction [24J and (b) 
the dipole-dipole interaction causes a blockade, due to which the probability 
for both atoms to be simultaneously excited (^44) reduces drastically. 



3.2 Non- identical atoms 

For the case of non-identical atoms, we assume different energy level separa- 
tion between ground and excited states of the two atoms and hence Ai 7^ A 2 . 
We also assume that the two decay constants 721 and 731 are close enough to 
be assumed equal which we denote by 7. The Rabi frequencies a 2 \ and a 3 i 
are also taken to be equal to a. Without a need to choose specific atomic 
species, we assume a situation of A 2 = Ai + 10, in units of 27, the population 
decay rat^j]. 

The steady state level populations for the case of non-identical atoms, as 
a function of the total detuning A = Ai + A 2 , are presented in figure [3J. In 
case of g = 0, it can be noticed that there are two resonances (dips) for the 



ground state population, as seen in the top curve of figure 3(a) These two 
resonances respectively occur when Ai = and A 2 = 0, and hence are on 
either side of the line-center A= 0. Corresponding to these two dips in the 
ground state population, a single resonance peak is seen in each of p 22 (at 
Ai = 0) and P33(at A 2 = 0). p 44 shows two very small peaks, again at ex- 
actly the same positions as the dips in p n . As the dipole coupling strength 
g is increased, side bands start to appear at +2a/ (\g\ 2 + 5 2 /4) + 2a 2 and 
—2a/ (\g\ 2 + <5 2 /4) — 2a 2 , where 5 = A x — A 2 . For higher values of g, the 
side bands are dominant only in p 22 and P33, corresponding to which there 
is hardly any activity in p 44 . The two-atom excited state p 44 , however looses 
the twin resonant peaks which were present in the absence of dipole - dipole 
interaction g, but develops a single peak at total detuning A = 0. This is in 
concurrence with the result of Varada and Agarwal p3|, indicating that as 
the strength of dipole-dipole interaction g increases, the probability for step- 



1 The two atoms need not be very different from each other. Instead of two different 
species, one could always use the same atoms, but two different transitions based on 
different polarizations and yet achieve the same effect. Although it is almost impossible 
to obtain two atoms such that their decay rates as well as dipole transition elements are 
exactly identical, one can choose a situation to begin with wherein the differences are 
insignificant compared to the energy level separation thus justifying our choice of various 
parameters above [25] . 
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wise excitation of both atoms decreases. At the same time, the probability 
for simultaneous excitation of both atoms through a two-photon absorption 
process increases due to the influence of g, as indicated by the central reso- 
nance peak of p 44 . Interestingly, the height of p 44 peak which is appearing 
only in the presence of g does not increase monotonically with an increase in 
g. There is an optimum value of g corresponding to a given value of the Rabi 
field strength a, at which p 44 attains its maximum value. Further increase in 
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g seems to result in a decrease in the peak height of p 44 , indicating a decrease 
in the probability for simultaneous absorption of two photons by the atoms, 
a feature which is present only for the case of non-identical atoms. The cor- 
responding peaks in p 22 and P33 at A = can be interpreted as due to decay 
of population from p 44 . The fluorescence from the system which is a linear 
combination of the level populations, namely p 22 + P33 + 2p 44 , shows quali- 
tative agreement with the experiment involving two terylene molecules [25J. 
We do observe an additional feature, namely the non-monotonic increase of 
the central peak as the dipole coupling strength is steadily increased. 

Figures 0] and |5] elucidate the above mentioned behaviour, specifically of 
the level populations, more clearly. In Figure HJ the population of the two- 
atom excited state p 44 is plotted for both identical (a) and non-identical (b) 
atoms for a value of the Rabi field strength a = 2 and for a range of values 
of the dipole coupling strength as indicated in the figure. There is a gradual 
decrease in the height of the central peak with increase in g for identical 
atoms. When it comes to the non-identical atoms, this behaviour is not re- 
peated and it is observed that the decrease in peak height with increase in 
dipole coupling strength is non monotonic which strongly suggests a compe- 
tition between a and g. To understand this behaviour for a continuous range 
of dipole-coupling strength, we have studied the on-resonance behaviour of 
p 44 (the behaviour of the central peak) as a function of the dipole coupling 
parameter for different values of the Rabi field strength g. The results are 
presented in Figure [5J where the on-resonance population of the two-atom 
excited state p 44 is plotted for different values of a as the dipole coupling 
strength is scanned from to 100. 

It is observed from figure 5(b) that for different values of a, there ex- 
ists a particular value of g for which p 44 attains its maximum and then it 
begins to decrease. This is not the case with identical atoms, wherein p 44 



monotonically decreases with increase in g as evident from figure 5(a) The 



non-monotonicity of the peak height is seen to occur only for non-identical 
atoms. 

The behaviour of the two-atom two-photon coherence Im p 14 is shown, 
both for identical as well as non-identical atoms in 6(a) and 6(b) respectively. 

The imaginary part of the atomic coherence p 44 [cf. figure Ej seems to 
exhibit non monotonic behaviour with respect to increase in g for both the 
cases. The height of the central peak is found to increase up to a value of 
g = 20 and a further increase in g(= 40) results in a decrease of the same. 
For the choice of detunings that are chosen in this figure, the side bands are 
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Figure 4: Level population P44 for the case of both (a) identical and (b) non- 
identical atoms for Rabi frequency a = 2 
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Figure 5: Level populations ^44 for (a) identical atoms \(b)\ non-identical 
atoms 



not well resolved. However, for large enough detunings, we did observe that 
these side bands are indeed well resolved. 
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In the results shown above, the atomic cooperative behaviour is incorpo- 
rated through the dipole-dipole coupling strength g. In the work of Varada 
and Agarwal [H], they have included both the cooperative decay parame- 
ter as well as the dipole-dipole coupling strength. We have carried out the 
same calculation and found that cooperative decay parameter acts like an 
additive feature to the dipole-dipole interaction strength. In all the cases 
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studied above, it is seen that the inclusion of the cooperative decay causes 
faster decay from level |4) resulting in enhancing the effect of g. But the fea- 
tures of non-monotonic behaviour of the level populations, coherences that 
are reported above remain unaffected due to the inclusion of the cooperative 
decay parameter. 

4 CONCLUSION 

We have studied the effect of dipole-dipole interaction on a pair of two-level 
atoms interacting with a laser field. The coupling of each of the atoms to the 
external radiation field, which is treated in the electric dipole approximation, 
and the strength of dipole interaction are incorporated phenomenologically 
by the parameters a and g respectively. The level populations and the in- 
duced atomic coherence are studied for a wide range of parameters both for 
identical as well as non identical atoms. The non-identical nature of the 
atoms is incorporated through different energy level separation while the 
other parameters such as the decay constants and coupling to the external 
field are assumed to be the same for both atoms, though in reality they may 
differ. Extensive studies are carried out for the above mentioned features of 
the atomic system for a wide range of parameters with particular emphasis 
on the competing effects of dipole coupling parameter and the Rabi field 
strength. In particular, the two photon resonance which appears in the pres- 
ence of the dipole coupling strength shows an increase upto certain values of 
the coupling strength and beyond the critical value, it decreases thus exhibit- 
ing a non-monotonic behaviour (only for the case of non-identical atoms). 
However, non monotonicity is seen in the coherence both for identical 
as well as non-identical atoms. The fluorescence signal shows qualitative 
agreement with the experimental results of ref [25] . 

Shaik Ahmed gratefully acknowledges the financial support from Univer- 
sity Grants Commission, Govt, of India. 
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